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ABSTRACT 



At the present time the primary method of obtaining solutions to 
non-linear differential equations is by means of the digital computer 
and numerical techniques. A method is here proposed to find an approxi- 
mate mathematical expression through the use of Laplace Transform tech- 
niques. Thus, the Laplace Transform concept is extended to the solution 
of non-linear differential equations. 
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CHAPTER 1 



1.0 Introduction 

Laplace Transforms have long been employed in the solution of 
linear differential equations. Baycura [ref. 1] proposed a method of 
obtaining an approximate solution to a non-linear differential equa- 
tion by using a linearized transform to obtain the approximation. In 
a later paper [ref. 2] it was proposed that an exact solution might be 
found by developing an expression for the Laplace Transform of a non- 
linear term. The result was that the solution was not exact, but very 
similar in form. 

These results are examined in detail here, and a general expres- 
sion for the non-linear transform is developed. Further, the concept 
is extended to non-homogeneous equations. 

1.1 Applications 

A method of obtaining a good approximation to the many non-linear 
differential equations which occur in all engineering fields would be 
of great utility to the engineer. Having a mathematical expression 
which closely approximates the solution would enable the engineer to 
analyze the system, and to optimize the design more rapidly, and 
without having to rely totally on digital computers. 
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CHAPTER 2 



2.0 Development of the Laplace Transform for Non-Linear Functions 
In order to use Laolace Transform techniques in obtaining an 
approximate solution to a non-linear differential equation, expres- 
sions for the transform must be derived. 



2.1 The Type of Equation Defined 

Consider a general non-linear differential equation of the form: 



where 

raised to the 



^ + a^ x^t) = f(t) (2.1.1) 

is the m^^ derivative of x(t); k indicates the derivative 
power; and x' (t) is x(t) raised to the n^^ power. 



2.2 The Method of Successive Integration by Parts 

To derive an expression for the transform, use is made of succes- 
sive integration by parts. Consider two time dependent functions, 
x(t) and y(t). The integral of their product is given by: 

/xy dt = X / y dt - X // y dt dt + x /// y dt dt dt - ••• 

+ (-l)"x^"^ / [/ y dt]" dt (2.2.1) 

where j]/ y dt]" is an n-fold iterated integral. This equation can 
now be used to find a series expression for the Laplace Transform of 
x(t) by lettinq y(t) = e“^^. Thus: 



-CO 

xe 

•*0 



-St 



dt 



X 



f<30 

e 



-St 



dt 



e“^^ dt dt + X 



.00 









-St 



dt dt dt - • 

( 2 . 2 . 2 ) 
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After integration the series becomes: 



fOO 

xe 



-St 



dt = - 




oo 





(2.2.3) 



In evaluating the integral at the upper limit the condition must be 
such that: 

Limit xe~^^ = 0 (2.2.4) 

t 00 



and similarly for all derivatives of x(t). Thus, evaluating the 
integral at the lower limit, the result is: 



Ce-^t . xM , iM , . ... 



(2.2.5) 



and since; 



X(s) 



,oo 

xe 



-st 



dt 



( 2 . 2 . 6 ) 



o 

by definition, an infinite series expression for the transform has 
been derived. 



2.3 Derivation of an Expression for the Transform of Derivatives 
To derive an expression for the transform of a derivative of 
x(t), the method just illustrated is used again, with x(t) in place 
of x(t). Thus: 



• « C “f • 

xe dt = X 



f-st 

e dt - X 



fe'^^ dt dt + X 



e"^^ dt dt dt - 



(2.3.1) 

Now, by adding and subtracting x(0), and then multiplying by s/s, the 
result is: 

x“^^ dt = s[^^ + - x(0) (2.3.2) 
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The term in brackets is just X(s), so: 

of[x(t)] = sX(s) - x(0) (2.3.3) 

a result which is well known. In a similar manner it is found that: 

c^CxCt)] = s^X(s) - sx(0) - x(0) (2.3.4) 

This method can be continued and it soon becomes apparent that the 

‘t* h 

expression for the transform of the m ^ derivative is: 

«f[x^'^^t)] = s^(s) - s"’'''x(0) - s"’'^x(0) - ••• - sx^"’"2^(0) 

- x^*^'^^(0) (2.3.5) 



2.4 Transform Expression for Integral Powers of a Function 

The functions considered thus far have all been linear; now the 
functions which cause the non-linearities will be considered. First, 



consider the function squared. 

-00 

x^e'^^ dt = |x(xe"^^)dt 



= X 



fOO . 

xe"^^ dt - X 



fl.-st 



dt dt + X 



.00 
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x'^^ dt dt dt 



(2.4.1) 



Now, let 



h 



r-st 

X dt 



(2.4.2) 



The integral is the same as Equation (2.2.6), so: 



II = x(0)X(s) 



Now let 



I 2 = - X 



'7 -St 

xe dt dt 



(2.4.3) 



(2.4.4) 



and if the inner integral is examined: 



xe'^^ dt = X 



e'^^ dt - X 



e"^^ dt dt + X 



-St 



-St *-st ’'-St 
X X _ X 

s " s2 ” 



dt dt dt - • • • 
(2.4.5) 
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This result is then put in Equation (2.4.4) to obtain: 






dt . 1, 



dt . 



roo . 

dt + 



(2.4.6) 



Each of the integrals can be evaluated by using the results of Section 
2.3, so: 






]_ 

S 3 



s^X(s) - sx(0) - x(0) 



sX(s) - x(0) 

(2.4.7) 

If this is expanded and like terms are collected, the result is: 

I 2 = (2.4.8) 

Continuing in this manner it is found that: 

^[x“(t)] = [^(0) +iM+iM+ ...jx(s) (2.4.9) 

But the term in brackets is just sX(s), so: 

^[xMt)] = sX^(s) (2.4.10) 

Thus, the transform of the square of a function has been derived. Now 
consider: 






^e'^^ dt = X 



x^e dt - X 



x=^e'^^ dt dt + 



(2.4.11) 



The first integral has already been found, so, from Equation (2.4.10): 

X x^e"^^ dt = x(0)sX^(s) (2.4.12) 

o 

The second integral is evaluated in the same manner as Equations 
(2.4.4) through (2.4.8). Thus: 

sXHs) (2.4.13) 



- X 



dt dt = 



The remaining integrals are evaluated in a like manner and the outcome 
is: 

r./^X ../nx ../nX ( 2 . 4 . 14 ) 



c([x = (t)] = + ••• 
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Once again, the term in brackets is recognized as being sX(s), so the 
expression can be written as: 

of[x'(t)] = s^xMs) (2.4.15) 

If higher powers of the function are examined it soon becomes apparent 
that: 

,iC[x"(t)] = s"-''x"(s) (2.4.16) 



2.5 Transform Expression for Powers of a Derivative 

Since many equations involve a derivative which is raised to a 
power, it will be useful to find an expression for the transform of 
such a function. First consider the square of the first derivative: 



(x)^e 



-st 



dt 



•xe"^^ dt - X 



xe'^^dt dt + • . 



(2.5.1) 



0 0 0 

The integral in the first term has been evaluated in Section 2.3, so 
this term becomes: 



xe 



■St 



dt = x(0)[sX(s) - x(0)] 



(2.5.2) 



0 

The second term is evaluated in the same manner as Equations (2.4.4) 
and (2.4.13). Therefore: 



- X 



fxe”^^ dt dt 



x(0) 



[sX(s) - x(0)] 



(2.5.3) 



0 

The remaining terms in the series are evaluated in the same way, so the 
transform becomes: 



= [sx(s) - x(0)][i(0) ..n 

= [SX(S) - x(0)] (s)[^ + ^+ •■■] (2.5.4) 
But the second bracketed term is equivalent to sX(s) - x(0), so the 
transform finally becomes: 



J^{[x(t)]^} = s[sX(s) - x(0)]^ 



(2.5.5) 
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It should be noticed that this last expression may be written as: 

(2-5.6) 

If the foregoing method is applied to the third power of the first 
derivative, the result is: 

?([x(t)]’} = > (2.5.7) 

If higher powers of the first derivative are examined, it can be 
shown that: 



x{[>i(t)]"} = s"-'/^[i(t)]}" 



(2.5.8) 

Now powers of the second derivative will be examined. First consider 
the second derivative squared: 



dt = X 



xe dt - X 



xe~^^ dt dt + 



(2.5.9) 



The integral in the first term is known to be s^X(s) - sx(0) - x(0). 

The second integral is treated in the same manner as Equations (2.4.4), 
(2.4.13), and (2.5.3) and becomes (l/s)[s^X(s) - sx(0) - x(0)]. Thus 
it is found that: 

.?/t;(t)]^} = [5(0) t ••■][s=x(s) - sx(0) - >i(0)] 

(2.5.10) 



The first term in brackets can be written as: 

s[hA + ‘ML+ ... 



= s' 



= s 



- sx(0) 

- sx(0) 



x(0 ) ^ x(0) ^ ‘x(0) ^ . 

s s 2 

(0) +MO) ... . x(o)‘ 

= ...j _ 

■ ■■■] ■ (2-5-n) 

This is recognized as being s[s'X(s) - sx(0) - x(0)]. Thus: 

T([x(t)]j = s[of[x(t)]]“ (2.5.12) 
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2.6 Summary of the Transforms 

If Equations (2.3.5), (2.4.16), (2.5.5), (2.5.7), and (2.5.12) are 
examined, it can be seen that a general expression can be written for 
the transform. Thus: 

( 2 . 6 . 1 ) 

where m = 0, 1, 2, 3, . . . ; and n = 1, 2, 3, ... 
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CHAPTER 3 



3.0 Application of the Transforms 

The use of the Laplace Transform expression just derived will now 
be demonstrated by obtaining solutions to several equations. 



3.1 Application to a Non-Linear, Homogeneous Equation 

Consider a homogeneous, non-linear differential equation of the 

form: 



e(t) + Ae(t) + BeMt) = 0; e(0) = V (3.1.1) 

Equation (2.6.1) is now used, and the transformed equation is: 

sE(s) - e(0) + AE(s) + BsE^s) = 0 (3.1.2) 

This expression is now rearranged to give: 

BsEMs) + (s + A)E(s) - V = 0 (3.1 .3) 

This equation is quadratic in E(s), so it may be solved for E(s) 
through the use of the familiar quadratic formula. Thus: 

E(s) = 



(s + A) ± \J(s + A)^ + 4BVs 
2Bi 



(3.1.4) 



Now this result is rewritten to put it into a more convenient and 
useful format. Thus: 



E(s) 



(s+A) 

^Bi~ 



■1 ± 



1 + 



4BVs 



Wat- 



(3.1.5) 



Now, if the positive radical is used and is expanded in a binomial 
series, the resulting expression for E(s) is: 

V BV^s . 2B^V^s^ 5B^V“s^ . 



E(s) 



BV^s , 2B^V^s^ 
TsW (s+Aj's" 



(3.1.6) 



s+^ ■ TsW ■ (s+A)5 ' Cs+A)’ ■ " 

Thus, an infinite series expression has been obtained for E(s). Now 
the inverse transform is to be found for each term of the series [ref. 3], 
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and the resulting expression is: 
.-At 



e(t) = Ve' 



1 - BVt 



1 M 
' ■ 2 



+ (BVt)2 



1 _ + AHM 



12 



- (BVt)^ 



f5 5At . AH^ AH^ 

6 ■ “¥~ “8 TW 



(3.1.7) 



An infinite series approximation has been found for e(t). However, for 
the original equation, an exact solution is known [ref. 4]: 



e(t) = Ve"^^ 



1 

' " A 



^e-At . 1 



-1 



(3.1.8) 



By expanding the term in brackets, the two solutions may be compared. 
Thus, after using a binomial expansion, the exact solution is: 



e(t) = Ve 



•At 



1 + 



BVf -At 



- 1 



2\/2 



B^V 



A2 



- 1 



3\/3 f 



B^V 






-At 



(3.1.9) 



Now, if each of the expotentials is expanded in a series, the result is 



e(t) = Ve 



■At 



1 - BVt 



At , AH^ AH" 



^ ■ 2 ' 6 






+ BHH^ 



1 - At + 



7AH^ 
12 



34.3 



A"t 



(3.1.10) 



Notice the similarity between this equation and the approximate solu- 
tion in Equation (3.1.7). Each sub-series in the approximate solution 
is truncated, since the second term of the s-domain expression yields 
two time-domain terms; the third term of the s-domain expression gives 
three time-domain terms; and so on. 

From an analysis of the parameters A, B, and V, it is apparent 
that BV < A, and further, the smaller BV is with respect to A the 
better the approximation will be. Table 1 shows a comparison of the 
exact and approximate solutions for two sets of B, V, and A. Figure 1 
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shows a graph of the exact and approximate solutions when A = 1.0 
B = 0.125, and V = 1 .0. 



TABLE 1 





A. 


A = 1 .0, B 


= 0.5, V = 1 .0 




Time 




Exact 


Approximate 


Error 


0.0 




1.000 


1 .000 


0.000 


0.5 




0.507 


0.473 


0.034 


1.0 




0.279 


0.257 


0.022 


1.5 




0.161 


0.174 


-0.013 


2.0 




0.094 


0.118 


-0.024 


2.5 




0.056 


0.057 


-0.001 


3.0 




0.034 


-0.010 


0.044 






A = 1.0, B = 


0.125, V = 1 .0 




0.0 




1.000 


1.000 


0.000 


0.5 




0.578 


0.577 


0.001 


1.0 




0.341 


0.344 


-0.003 


1.5 




0.203 


0.214 


-0.011 


2.0 




0.122 


0.139 


-0.017 


2.5 




0.074 


0.095 


-0.021 


3.0 




0.044 


0.067 


-0.023 


3.5 




0.027 


0.048 


-0.021 


4.0 




0.016 


0.035 


-0.019 


4.5 




0.010 


0.026 


-0.016 


5.0 




0.006 


0.019 


-0.013 
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i o 

f • 

1 
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3.2 Application to a Non-Linear, Non-Homogeneous Equation 
Consider an equation of the form: 

Mv + kv^ = Mg (3.2.1 ) 

This is the equation of motion of a falling object with the retarding 
effect of air resistance considered. It is assumed that the retarding 
force is proportional to the velocity squared. If the equation is 
written as: 

V + hv^ - g = 0 (3.2.2) 

where 

v(0) = 0 (3.2.3) 

and 



h = k/M, (3.2.4) 

then an approximate solution to the equation can be found by applying 
Equation (2.6.1) as before. Thus: 

sV(s) - v(0) + hsVMs) - 1= 0 (3.2.5) 

This expression is rearranged to get: 

hs^V^(s) + s^V(s) - g = 0 (3.2.6) 

Once again the result is a quadratic equation in V(s), and it is 
solved for V(s) by using the quadratic formula. Thus: 

V(s) = s^ ± ^s“ + 4ghs^J (3.2.7) 



This equation is now rearranged to get the same form as that used in 
Equation (3.1.5): 



V(s) 



1 



1 ± 



^ 4^1 

S^ 



(3.2.8) 



If the term in parentheses is expanded in a binomial series and the 
positive sign in front of the radical is used, the result, after 
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simplification, is: 



V(s) ••• (3.2.9) 

and the time-domain solution is obtained by finding the inverse trans- 
form of each term of the series. Thus: 

:: ., 4 . q^ht^ ^ 2q^h^t® 5q‘‘hH^ 14q®h‘*t® tr, o ir,\ 

v(t) = qt - • • • (3.2.10) 

/ 

So an approximate solution has been found. The exact solution is 
known [ref. 5]: 



v(t) = 





^ _ g-2(gh)^t 


F 

V / 


T + e‘2(gh)^t 



(3.2.11) 



where (q/h)'^ is the terminal velocity of the falling object. If the 
numerator and denominator of this expression are multiplied by e^^ ' 
then: 



v(t) = 



^ tanh (qh)^t 



(3.2.12) 



The hyperbolic tangent can be written as an infinite series, and in 



this form the exact solution is 

r— 3 



V(t) = - Istipi 2(2!^. ...j (3. 2. 13) 

To compare the two solutions. Equation (3.2.10) is rewritten as: 
v(t) = 





g 


% 


(gh)"^ - 


pl(9t!)“t’ . 


ni 


2(qh)'t^ f 5 ]l7(qh)H^ ... 




h 






8 

k. J 


^T5 3Tt 



(3.2.14) 



The series are similar in form, but each term of the approximate solu- 
tion is smaller than the corresponding term in the exact solution. In 
order to get a numerical comparison, consider the falling object to be 
a spherical rock with density (p) of 2600 kilograms/meter^. The 
following parameters must be defined and evaluated [ref. 6]: 

M = pV 

k = 0.5 p^ Cj A 
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where M is the mass of the object, V is its volume, and A is its 
surface area. The other constants are: 

= 1.293 kilograms/meter^ (Actually, is a function of air 

pressure and temperature, but it 
will be assumed to be very nearly 
constant for this analysis.) 

Cj = 0.5, the drag coefficient for a sphere 

A = ttD^, the surface area of a sphere 
ttD^ 

V = -g— , the volume of a sphere 

g = 9.8 meters/second^ , the acceleration due to gravity 
Thus, if the sphere has a diameter of 0.1 meter, the exact solution is: 

v(t) = 36.4 tanh 0.27t (3.2.15) 

and the approximate solution becomes 



v(t) s 36.4 




(0.27t)^ . (0.27t)5 
6 60 



(0.27t)’ . (0.27t)® 

1008 25920 ■ ■ * 

(3.2.16) 



The comparison of these two expressions is shown in Table II and a graph 
in Figure 2. In the graph, notice that the exact solution asymptoti- 
cally approaches the terminal velocity 36.4 meters/second, while the 
approximate solution overshoots that value, and after about seven 
seconds decreases monotonically. This is due to the fact that the 
highest power term in the series was preceded by a negative sign, and 
as t became large this term overpowered the others. The first term 
omitted after the series was truncated was positive, so if it were 
included, the graph would end by increasing monotonically. The hump 
is caused by the lower degree terms being predominant when t is small. 



19 



TABLE II 



Time 


Exact 


Approximate 


Error 


0.0 


0.000 


0.000 


0.000 


0.5 


4.884 


4.885 


-0.001 


1.0 


9.596 


9.682 


-0.086 


1.5 


11.839 


12.021 


-0.182 


2.0 


17.945 


18.679 


-0.734 


2.5 


21.413 


22.731 


-1.318 


3.0 


24.373 


26.403 


-2.030 


3.5 


26.845 


29.649 


-2.804 


4.0 


28.872 


32.440 


-3.568 


4.5 


30.510 


34.757 


-4.247 


5.0 


31.816 


36.598 


-4.782 


5.5 


32.847 


37.968 


-5.121 


6.0 


33.656 


38.879 


-5.223 


6.5 


34.287 


39.337 


-5.050 


7.0 


34.776 


39.328 


-4.552 


7.5 


35.153 


38.788 


-3.635 


8.0 


35.444 


37.574 


-2.130 


8.5 


35.668 


35.399 


0.269 


9.0 


35.840 


31.750 


4.090 


9.5 


35.972 


25.768 


10.204 


10.0 


36.073 


16.077 


19.996 
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Velocity vs Time 

Comparison of the exact and approximate solutions for 
the equation: Mv(t) + kv^(t) = Mg 
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Figure 2 seconds 



3.3 Application to an Equation with a Third Degree Non-Linearity 
Consider now a non-linear differential equation of the form: 
x(t) + w^x(t) + hxMt) = 0, x(0) = A, x(0) = 0 (3.3.1) 

After using the transform expression of Equation (2.6.1), the trans- 
formed equation is: 



hs"X3(s) + (s"+to2)X(s) - As = 0 (3.3.2) 

This is in the so-called "normal" form of a cubic equation, in that 
there is no term involving X^(s) present. Thus there are solutions 
X (s), X (s), and X (s); where: 

12 3 

X = A + B, and X , X = - %(A+B) ± ^ (A-B) (3.3.3) 

1 2 3 

and where: 



A,B 




(s2+(4)3 

27h3s6 



(3.3.4) 



It can be seen that expanding this expression and determining which 
root to use is a formidable task. However, it is known from Equation 

(2.2.5) that: 

X(5) ... (3.3.5) 

and if the initial conditions from Equation (3.3.1) are used and if it 
is assumed that all higher derivatives can be neglected at t = 0, then 
Equation (3.3.2) may be written as: 

(3.3.6) 

Thus, the degree of the equation has been reduced by one and the solu- 
tion can be obtained in the same manner as in the previous two sections 

4hA^s^ 



hAsX^(s) + (s"+a)J)X(s) - As = 0 



X(s) = 

2hAS 



■1 ±1 + 



(s2+o)J)== 






(3.3.7) 



Then, expanding the term in brackets as in previous solutions, the 
result is: 

As hA^s^ , 2h^A®s® , in n q\ 

" s2+^ “ (s2+(j 2)3 + (sHoiJ)^ (3.3.8) 
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Now the inverse transform is found for each term of the series. For 
the first term, then: 



-1 



As 






= A cos CO t 
0 



(3.3.9) 



The inverse transforms for the succeeding terms are more difficult to 
evaluate, as no expression of the form s'^/(s^+Wq)’^ can be found in the 
tables of Laplace Transforms available. However, one handbook [ref. 7] 



(3.3.10) 



does have the transform: 

Further, the table [ref. 8] has the relation that: 

s"g(s) = f^^ht) + f^'^^(0)s'^”^''^, n = 1, 2, 3, ••• (3.3.11) 

k=o 

f k) 

if f' '(0) =0 for k = 0, 1 , . . . , n-1 . Thus, these two relations can 
be combined to yield the desired inverse transforms. So, for the 



second term of the series: 
s 



g(s) = 



(S^+U)^)^ 



(3.3.12) 



and its inverse transform is given by: 



f(t) 



2^r(3) 



(3.3.13) 



to. 



Then, the Bessel Function is given by [ref. 9]: 

rsin cDot 



J-2.(o)-t) = 
2 0 



2 1 


n 


[iTco^tJ 


1 [ 



CO t 
0 



- cos co^t 
0 



(3.3.14) 



Upon substituting this expression into Equation (3.3.13), the resultant 
inverse transform of the second term is, after applying Equation 
(3.3.11): 



f(t) = (3 10 t sin 10 t +co^t^ cos CO t) 

' 8co? o 0 0 0 

0 



(3.3.15) 
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This method is then applied to the other terms of the series and the 
resulting expression is: 

hA^ 



x(t) = A cos y - (3y sin y + y^ cos y) + 

5y2 . . 5''^ 



h^AMS 

Ipy sin y 



■ ^ sin y + cos y 



(3.3.16) 



where, for convenience, y = lo^t. Another approximate solution is given 
by Cunningham [ref. 10], using the Ritz method. The result is: 



x(t) = A cos'\jo^ + 0.75hA^ t 



(3.3.17) 



Using the arbitrary values = 1 > h = 1 , and A = 1, a comparison of 
these two solutions is shown in Table III and a graph in Figure 3. 

With the arbitrary values selected, the new method yields good results, 
when compared to the Ritz method, up to about one second. Examination 
of the series shows that better results could be obtained over a longer 
period of time if the quantity A were made less than unity. 

In the method originally proposed by Baycura [ref. 1], the trans- 
formed equation was the same as that which would result if the relation 
X(s) equals A/s were used again in Equation (3.3.2): 

X(s) = (3.3.18) 

For this expression, then the inverse transform is: 

x(t) = cosA/2 t (3.3.19) 

assuming that A = h = = 1 . At t = 2, x(t) = -0.951, and at t = 1 , 

x(t) = 0.156. Thus, if desired, the equation in s can be reduced to a 
linear expression, and this will yield a solution which is not in the 
form of an infinite series. 
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Displacement vs Time 

Comparison of two approximations to the 

solution of: x(t) + ^x(t) + hx^t) » 0 



O 

CO 
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Figure 3 seconds 



TABLE III 



Time 


Ritz 

Method 


Laplace 

Method 


Error 


0.0 


1 .000 


1 .000 


0.000 


0.2 


0.965 


0.960 


0.005 


0.4 


0.863 


0.844 


0.019 


0.6 


0.701 


0.658 


0.043 


0.8 


0.490 


0.416 


0.074 


1.0 


0.245 


0.134 


0.111 


1.2 


-0.017 


-0.173 


0.156 


1.4 


-0.278 


-0.485 


0.207 


1.6 


-0.519 


-0.786 


0.267 


1.8 


-0.724 


-1.060 


0.336 


2.0 


-0.880 


-1.290 


0.410 



3.4 An Equation with a Derivative Raised to a Power 
Consider now an equation of the form [ref. 11]: 
x(t) + [x(t)]^ + x(t) + k = 0, x(0) = 1, x(0) = 0 
The transformed equation is found as before, thus: 



s^X2(s) - (s2-l)X(s) +1=0 

This is then solved for X(s) as before, and the result is: 



X(s) 



- (s^-1) 

" 2s3 



1 ± 



1 - 



4ks^ 



(3.4.1) 



(3.4.2) 



(3.4.3) 
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The term in brackets is then expanded as a binomial series, but this 
time the negative sign in front of the radical is chosen to avoid a 
delta function in the solution, so that: 



_ k k^s^ . 2k^s® . 

+ (s2-l)3 (s2-1)5 + •• 



(3.4.4) 



The inverse transform of the first term is k sinh t, but for the suc- 
ceeding terms the same difficulty is encountered as in Section 3.3. 
However, the table of transforms [ref. 12] lists: 






(S2-1)' 



I 3 



»-2(t) 



2"‘^r(n) 



(3.4.5) 



where I^_^(t) is the modified Bessel function [ref. 13] given by: 
2 



^N(t) 



= I 



N-^2k 



k=o 



(t/2) 
k!r(N-Hk+l ) 



For the second term of the series, n = 3, N = 3/2, and: 

= J mi? 



1 1 
2 



. (t/2)2 (t/2)f (t/2)' 

■Tdr 



and to evaluate the Gamma function [ref. 14]: 



T>/ ii\ 1.3.5.7...( 2m- 1 ) I — 

r(m+J$) = fjp -^Jrr 



(3.4.6) 



(3.4.7) 



(3.4.8) 



Thus: 



1 1 



I (t) = ^(t/2)^ ^ 8(t/2)^ ^ 16(t/2) 

2" 3-\fW 1 Bi/tT 1 05-^ 

This result is substituted in Equation (3.4.5) to get: 

cT' 



1 


s 


1 


r4t“ . t®, t® , 

- - _ -j- M ' m m-i *r • • • 




U*-l)^ 


3? 


^ T? TT5F 



(3.4.9) 



(3.4.10) 



Examination of this expression shows that the system is unstable, since 
the terms are always positive and the series has no limit, increasing 
without bound as t becomes large. 
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CHAPTER 4 



4.0 Conclusions 

It has been established that a method for obtaining an approximate 
solution to a non-linear differential equation has been derived, and 
that the method uses techniques and characteristics inherent to the 
Laplace Transform. Further, it has been shown that sufficient accuracy 
exists in the approximation to make it useful in many engineering 
appl i cations. 

As in any approximation which uses a truncated infinite series to 
obtain a numerical answer to a problem, extreme care must be exercised 
to ensure that the approximation is accurate within a reasonable degree. 
This means that the parameters must be examined closely to avoid 
exceeding the capability of the approximation. It also means that the 
approximation may indeed only be accurate over a small part of the 
entire solution. 

The results included here are the outcome of what really amounts to 
preliminary investigation only. One question that remains is that, 
since the solutions obtained in Sections 3.1 and 3.2 are so similar in 
form to the exact solutions, and since the solution seems exact in the 
Laplace expression, why isn't it, in fact, exact? After all, one major 
attraction of the Laplace Transform in solving linear equations is that 
the solution may be obtained by manipulating the transformed equation 
algebraically, and that is how the solutions were arrived at here. Thus, 
it seems that further research must be done in the area of this non- 
linear transform. 
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